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Non-Markovian evolution of an open quantum system can be 'unraveled' into pure state trajec- 
tories generated by a non-Markovian stochastic (diffusive) Schrodinger equation, as introduced by 
Diosi, Gisin, and Strunz. Recently we have shown that such equations can be derived using the 
modal (hidden variable) interpretation of quantum mechanics. In this paper we generalize this the- 
ory to treat jump-like unravelings. To illustrate the jump-like behavior we consider a simple system: 
A classically driven (at Rabi frequency Q) two-level atom coupled linearly to a three mode optical 
bath, with a central frequency equal to the frequency of the atom, lvq, and the two side bands have 
frequencies ujq ± Q. In the large Q limit we observed that the jump-like behavior is similar to that 
observed in this system with a Markovian (broad band) bath. This is expected as in the Markovian 
limit the fluorescence spectrum for a strongly driven two level atom takes the form of a MoUow 
triplet. However the length of time for which the Markovian-like behaviour persists depends upon 
which jump-like unraveling is used. 
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I. INTRODUCTION 

In the past few years non-Markovian stochastic 
Schrodinger equations (SSEs) for diffusive unravelings 
have received considerable attention [1, 2, 3, 4, 5, 6, 7, 
8, 9, 10, 11, 12]. A diffusive non-Markovian stochastic 
Schrodinger equation is a non-hnear, non-Markovian, and 
stochastic evolution equation for a quantum state. We 
usually denote this quantum state by \ijj[zkit)}(t)) where 
the subscripts {zfc(t)} implies that \4'{zk{t)}{t)) is condi- 
tioned on the set of continuous time dependent random 
variables {zk(t)} described by a probability distribution 
P{{zk},t). This set is constrained such that 

PrcdW = E[\^{Mt)}{t)){^[,,^t)}m, (1.1) 

where E[...] denotes an ensemble average over the set 
{zk{t)}. Here prcd(i) refers to the reduced state of an 
open quantum system. An open quantum system is a 
combined system consisting of a "system of interest" and 
an environment or bath. The reduced state is defined by 

Prcd(t) =Trenv[|*(t))(*(t)|]. (1.2) 

Here \'^{t)) is the solution of the Schrodinger equation 
for the combined system (generally an entangled state) 
and Tronv[...] refers to a partial trace performed over the 
environments degrees of freedom. 

In this paper we extend the known unravelings for non- 
Markovian SSEs to include two jump-like unravelings. 
This is motivated by the fact that in the Markovian limit 
the simplest of all Markovian SSEs is jump-like in na- 
ture. This corresponding to the direct detection quan- 
tum trajectory [13]. However, as shown in Ref. [8] in 
the non-Markovian case SSEs are not quantum trajecto- 
ries. That is, they are not evolution equations for the 
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state of the system conditioned on the outcomes of some 
continuous-in-timc measurements of the bath. This is 
because a non-Markovian system remembers what has 
happen in the past (whether a measurement has occurred 
or not) and a continuously measured system will have a 
drastically different average evolution, p{t), from a non- 
measured system Prcd(i). 

So what is the interpretation of non-Markovian SSEs? 
In Ref. [11] we came to the conclusion that the three 
current non-Markovian SSEs, the coherent [2, 5, 8], 
quadrature [7, 8] and the position [12] only have a non- 
trivial (non-numerical) interpretation when we consider 
the modal [14, 15, 16, 17, 18, 19, 20, 21] (or Beh's be- 
able [22]) interpretation of quantum mechanics. They 
are evolution equations for the system part of the prop- 
erty state of the universe when certain bath observablcs 
are given an objective reality (become hidden variables). 
The values of these bath hidden variables, at time t, are 
given by the set {zk{t)}. The different unraveHngs corre- 
spond to a different choice of which bath observables are 
to be objective real. For example the position unraveling 
arises when the hidden variables for the bath correspond 
to the positions of each oscillator making up the bath 
[11, 12]. That is, the set {zk{t)} becomes the set {xk{t)} 
corresponding to the values of the positions of each oscil- 
lator of the bath. For this case the position trajectories, 
{(itXfc(t)}, actually correspond to Bohmian trajectories 
[23]. 

Note this is not the only interpretation of non- 
Markovian SSEs; in Rcfs. [7] and [24] Bassi and Ghirardi 
interpreted these equations as a new dynamical reduc- 
tion modal, dynamical reduction with non-white Gaus- 
sian noise. In this paper we are not going to argue against 
this interpretation except to say that by accepting the 
above modal interpretation we can easily generalize the 
theory for diffusive non-Markovian SSEs to include jump- 
like "unravelings." 

A difference from the diffusive case is that it does not 
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appear possible to derive an explicit expression for the 
jump-like non-Markovian SSEs. This is because these 
non-Markovian SSEs are conditioned on a discrete set 
{z„^ (t)} and as such it is impossible to define derivatives 
with respect to this set. Furthermore a Girsanov trans- 
formation [2] can not be defined. In fact it turns out that 
to derive numerical jump-like non-Markovian unravelings 
is effectively the same as solving the Schrodinger equation 
for the combined system. Thus any numerical advantage 
of using non-Markovian SSE to determine Prod(0 [see Eq. 
(1.1)] is lost for these unravelings. However, the point of 
this paper is not to derive more effective ways of finding 
Prcd(i) but rather to illustrate that by using the modal 
interpretation we can derive the numerical equivalent of 
the solution to a non-Markovian SSE for jump-like un- 
ravelings. 

To illustrate a non-Markovian unraveling that ex- 
hibit jump-like behavior we consider the two unravelings, 
spectral-mode and temporal-mode. We apply our theory 
to a simple system: A classically driven two level atom 
(TLA), at Rabi frequency fl, coupled linearly to a three- 
mode optical bath; with a central frequency equal to the 
frequency of the atom, uiq, and the two side bands have 
frequencies ujq ± fl. In the large fl limit we observed 
that the jump-like behavior is similar to that observed in 
Markovian systems. This is expected as in the Markovian 
limit a strongly driven TLA has a fluorescent spectrum 
which takes the form of a MoUow triplet [25] with peaks 
at these frequencies. 

The structure of this paper is as follows: In the next 
section we present an outline of the general physical 
model to which the theory of non-Markovian SSE is ap- 
plicable. In Sec. Ill we present a review of the modal 
interpretation of quantum mechanics and how it applies 
to non-Markovian SSEs. In Sees. IV and V we de- 
rive two jump-like unravelings, these being the spectral- 
mode and the temporal-mode unraveling. In Sec VI we 
briefly outline how Jack, Collet and Walls [26, 27, 28] 
non-Markovian quantum trajectory theory differs to our 
theory. Note other jump-like non-Markovian unravelings 
exists [29, 30], but these are only numerical methods for 
solving pred{t), and as such will not be considered in this 
paper. Finally in Sec. VII we conclude this paper. 



II. THE TLA AND THE UNDERLYING BATH 
DYNAMICS 

The aim of this section is to outline the underly- 
ing model we will be using to generate jump-like non- 
Markovian unravelings. The standard model is to assume 
a preferred factorization of the universe into a system and 
environment, label by the Hilbert spaces Tisys and Tionv 
respectively. The total Hamiltonian for the universe is 
then given by 

Hsys{t)<»i + i®Honv + V. (2.1) 



Here -ffsys, i?env and V represent the Hamiltonians for 
the system, environment and any interaction taking place 
between the system and environment respectively. In this 
paper we assume the system is a two level atom (TLA) 
described by the free system Hamiltonian 

Hsysit) = i?.yso + H (t) = + H (t) , (2.2) 

where ujq is the atomic transition frequency and is the 
Pauli spin operator 

a, = \e){e\-\b){b\, (2.3) 

where |e) and \b) represent the excited and ground state 
of the atom. The extra term H{t) in Eq. (2.2) repre- 
sents any extra system evolution, for example a classical 
driving process. 

The environment, as in all optical situations, is mod- 
elled by a collection of one-dimensional harmonic oscilla- 
tors. The free Hamiltonian for this type of environment 
is 

K 

Hcnv = ^ hujkalcLk, (2.4) 

k 

where k defines the total number of modes, while LUk and 
cik {a\) are the frequency and the annihilation (creation) 
operator of the k^^ mode of the environment. 

We assume that the interaction between the system 
and environment is consistent with the rotating wave ap- 
proximation. That is, we assume V is linear in the envi- 
ronment amplitudes, and has the form 

K 

V = i'h'Y2{glcra\- gko-^ CLk), (2.5) 
k 

where gk is the coupling strength of the fc**^ mode to the 
system. 

For calculational purposes we define an interaction 
frame such that the Hamiltonians -ffsyso [which is defined 
implicity in Eq. (2.2)] and i?onv is removed. The unitary 
evolution operator for this transformations is 

Uo{t, 0) = e-*(^-=o®i+i®^-v)(t-o)/n^ (3.6) 

Thus the combined state in the interaction frame is de- 
fined as 

|vl/(i)) = ;7t(f,0)|*(t)sch), (2.7) 
and an arbitrary operator A becomes 

A,^t{t) = ul{t,Q)AU^{t,Q). (2.8) 
This allows us to write the Schrodinger equation as 

dt\^{t)) = -'Mntit) + Vint(i)]l*W), (2.9) 
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where Hint{t) refers to H{t) in the interaction frame and 
Vint{t) is defined as 

Vlnt(t) = i?i^[<7g^e'"''*4 - atg,e-^"^*afc], (2.10) 

k 

where flk = ujk — uJo- 

In aU the examples we present in this paper we take the 
extra system Hamiltonian, H{t), to be classical driving 
of the TLA. Under the dipole and rotating wave approx- 
imation the extra system Hamiltonian in Eq. (2.2) is 

H{t) = ^[a'' cxp(-iuJct) + a exp{iujct)], (2.11) 

where ft = — 2d • Eo(ro)/?i is the Rabi frequency and 
Eo(ro) is the amplitude of the classical field at the posi- 
tions of the atom and d is dipole transition matrix. For 
simplicity we tune the frequency of the classical driving 
field, LUc, to the atomic transition frequency. Thus when 
moving to the interaction frame Eq. (2.11) becomes 

Hint = — (a^ + a) = —a,, (2.12) 
which is time independent. 

III. MODAL INTERPRETATION OF 
QUANTUM MECHANICS 

A. General modal dynamics 

In this section we give a brief overview of the modal in- 
terpretation of quantum mechanics; for a more detailed 
description see Rcfs. [17, 18, 19, 20, 21, 22]. The ba- 
sic idea of this view of quantum mechanics is that cer- 
tain observables have an objective reality independent of 
measurement. This is in contrast to the orthodox inter- 
pretation where reality is undefined prior to measurement 
and it is the act of observation that defines the reality of 
the observable. 

In quantum mechanics it is convenient to define an 
observable as 

Z = {(z„,7r„)}, (3.1) 

where the set of values {zn} correspond to the possible 
values the observable can have and {•S'„} is the projec- 
tive measure for the observable. When these values cor- 
respond to numbers we can write an observable as an 
operator 

Z = ^Zn'iTn- (3.2) 

n 

The projectors are orthogonal and form a decomposition 
of unity: 

^7r„ = i. (3.3) 

n 



In quantum mechanics because of the non- 
commutative nature of observables, not all observables 
can be written in terms of the same projective measure 
{tTji}. In the modal theory one postulates that only one 
set of projectors has an objective reality, this being called 
the preferred projective measure. Once this measure is 
specified it uniquely defines a group of observables which 
are objectively real, this being the group defined to have 
elements given by Eq. (3.1). In this paper we will refer 
to this group as the properties of the system. That 
is, the group of preferred observables (observables with 
objective reality) are labelled properties of the system. 

The biggest criticism against the modal theory is that 
many choices of the preferred measure are viable [17]. 
This has lead to many variants of the modal interpre- 
tation of quantum mechanics [15, 16, 17, 18, 19, 20, 
21, 22, 23]. Some have tried to address the problem 
of choice, either by accepting it (the beable variant) 
[17, 19, 22], selecting the position projective measure 
as preferred (Bohmian mechanics [23]), or by applying 
an algorithm for determining the preferred projective 
measure based on the total wavefunction of the universe 
[15, 16, 17, 18, 21]. The algorithmic approaches in our 
opinion still face choice as one has to choose the preferred 
factorization of the universe. Furthermore we have shown 
in Rcf. [20] that we can extend the modal theory to in- 
clude positive operator measures, thereby increasing the 
amount of possible choices. Here, however, we will only 
consider preferred projective measures and take the view 
that choice is fundamental. Depending on the physical 
situation we wish to describe we will choose the appro- 
priate measure {7r„}. 

Once we have chosen the preferred projective measure, 
to explain why property Z has the value Zn at time t we 
can introduce an extra quantum state, the property state. 
It is defined as 

\^.„it){t)}=n^\^{t))/VN, (3.4) 

where iV is a normalization constant and z„{t) is the 
value of Z at time t. The property state at time t is deter- 
mined by both |5'(t)) and a stochastic evolution (jumps 
between different n) . It is interpreted as the actual state 
of the universe, and by the eigenstate-eigenvalue link it 
selects the present value Zn for property Z from the possi- 
ble values {zn}. The stochastic dynamics (rates at which 
it jumps between difi'erent n) is determined by \^{t)) and 
as such |4'(t)) in this interpretation is called the guiding 
state. 

The modal dynamics (the stochastic evolution of the 
property state) is found using the method originally pro- 
posed by Bell [22] and gcnerahzcd in Rcfs. [18, 19, 20]. 
We start by defining Pr(z„, t) as the probability that the 
property will have the value z„ at time t. Assuming a 
Markovian process, by which we mean that the probabil- 
ity of the property being z^ at time t + dt only depends 
on the value at time t, we can write a master equation 
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for Pr(z„, t) as 
dtPr(z„,t) = ^[T,„„(t)Pr(z„„i)-r„„Pr(z„,t)], (3.5) 

m 

where Tnm (for ri ^ to) arc transition rates. For n = m, 
Tnn (which is negative) is a measure of the rate at which 
value z„ foses probabihty. 

Defining a probabihty current Jnm{t) as 

Jnm{t) = Tnm{t)Pr{Zm,t) - TmnP^^Zn, t) , (3.6) 

results in Jnm{i) — —Jmn{t) and allows us to rewrite the 
probability master equation as 



dtPr(z„, i) = ^ Jnm{t)- 



(3.7) 



Given Jnm{i) and Pr(z„, t), there are many possible tran- 
sition rates satisfying Eq. (3.7). One solution, chosen by 
Bell [22] is as follows. 
For Jnm{t) < 0, 

Tnm{t) = 0, (3.8) 

Tmn{t) = -J„™(i)/Pr(^n,i), (3.9) 

and for Jnm{t) > 

Tnmit) = J,™(t)/Pr(z„„t), (3.10) 
T,„„(t) = 0. (3.11) 

Thus once we have the probability current it is possible 
to calculate the transition matrix Tnmit) which in turn 
allows us to calculate (via using a random number gen- 
erator) a trajectory for the value of the property Z. 

To find Jnm{t)i as in the orthodox theory, we have to 
postulate a fundamental rule for probability, the Born 
rule 



Pr(z„,t) = (1'(t)|^„|*(0). 



(3.12) 



With this equation and Eqs. (3.7) and (2.9) a possible 
solution for J„,„(i) is 

Jnm{t) ^ 2Im{(*(i)|7r„[i?i„t(i)+T/i„t(t)] 

y.TTm\^{t))}/h. (3.13) 

Note the ensemble set of trajectories we obtained for the 
value of the property Z is only one of the infinitely many 
possible sets of ensemble trajectories. Others can be 
found by either adding an extra term T,j^(<) to r„„j(t), 
where T^^it) is constrained only by 

T°m{t)Pm{t) - T^n{t)Pn{t) - 0, (3.14) 

or by adding any current J„°„(i) to Jnmit) which satisfies 

m 

For the purposes of this paper we only consider Bell solu- 
tion [not containing the extra T,j^(t) and J„m(t) terms]. 
For a discussion of these solutions see Refs. [31] or [32]. 



B. Application to non-Markovian SSEs 

In Ref. [11] we showed that the theory of the non- 
Markovian SSEs emerge from modal dynamics when we 
assume a preferred projective measure of the form 



{■kn = Isys CS" 7rn„,„}j 



(3.16) 



where 7r„^j^^ is a projector define solely to operate in the 
Hilbert space of the environment. (Note the coherent 
non-Markovian SSE [1, 2, 3, 4, 5, 8, 11] arises when we 
use a preferred POM rather than a projective measure; 
see Ref. [11]). This means the observables which have 
definite values are of the form 



Z = {zn, i 



sys 



,}■ 



(3.17) 



That is, the bath is given property status, while the sys- 
tem is treated as a purely quantum system, which never- 
theless influences the bath values via the coupling Hamil- 
tonian [Eq. (2.10)]. 

Diffusive non-Markovian arise when the projective 
measure for the environment {t^u^^SS is chosen to be 
an infinitesimal projective density measure {T^[z)e^^dz^. 
That is, properties are given by 



{z, isys (g) 7r(z)cnvd2}, 



(3.18) 



where Z now has a continuous set of possible values 
{z}. For example the position unraveling occurs when 
we choose the following infinitesimal projective density 
measure 

{7r(a;i)c?xi ® ... ® Tt{xk)dxk ® ... (8) 7r(a;K)c?x„}, (3.19) 

where 7r(a;fc) = and \xk) is defined to be the 

eigenstate of the operator x/^ = (a^ + al.)/\/2. That 
is, the observables correspond to the positions of the en- 
vironment's K harmonic oscillators (and any function of 
them by the principle of property compositions [18]) are 
given property status. With this choice of preferred pro- 
jective measure the property state of the universe be- 



comes 



l^{..(t)}W> = \{xk{t)}) <E> |V'{..(t)}(i)), (3.20) 

where li^^Xkyit)) is called the conditional system state. 
Here we see that the property state is divided into two 
parts, an environment state and system state. 

In Ref. [11] we showed that when using Bell's solu- 
tion for the transition parameters we can derive a set of 
K trajectories for the values of the k position properties, 
which we denote as {xk{t)} (these actually turn out to 
be Bohmian trajectories). Using these trajectories we can 
derive a differential equation for |V'{xfc(t)}(0): which de- 
scribes the evolution of the system part of the property 
state of the universe, the environment part is given by 
|{a;fe(t)}). This system-state differential equation turns 
out to be the non-Markovian SSE for the position unrav- 
eling. For the complete derivation of this non-Markovian 
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SSE and the other two diffusive non-Markovian SSEs see 
Refs. [2, 8, 11, 12]. 

Here we are going to use the above reasoning to moti- 
vate how one can go about finding numerically the equiv- 
alent to the conditioned system state for jump-like un- 
ravelings. To do this we simply choose discrete rank one 
projectors for the environment. That is, the preferred 
projective measure is 



{7r„} = {isys® |{^«J)({2nJ|}, 



(3.21) 



where K^nJ) = |z„i) (g) ... (g) \zn^) <E) ... (E) \znj and n = 
{ni, 712, With this preferred projective measure 
the property state of the universe [sec Eq. (3.4)] becomes 

^ \Mt)}) ® l^{.„,(t)}W), (3.22) 
where the conditioned state is 

= {{znAtmit))/^- (3.23) 

The evolution of this state through time is equivalent to 
a non-Markovian SSE for jump-like unravelings. 

To summarize, the procedure to developed numerical 
trajectories for both |i/'{z„j^(t)}(0) ^^^^ {zn^it)} is to first 
calculate the guiding state from the Schrodinger equa- 
tion [see Eq. (2.9)], then using Bells solution we can 
determine both Jnm (t) and Tnm (t) for the appropriately 
chosen preferred projective measure. Once these are cal- 
culated we use a random number generator to simulate si- 
multaneously a typical trajectory for {znf.{t)} and hence 
I^WJt)}(0>. 



IV. THE SPECTRAL MODE UNRAVELING 

The first unraveling we consider is the simplest of all; 
it corresponds to the bath hidden variables being the 
photon number of each bath spectral mode. That is. the 
values of the bath hidden variables are denoted by {n^} 
and only take on discrete values. The preferred projective 
measure for this unraveling is 



(4.1) 



where \{nk}) = \ni) 



Ink)®. 



[n^) and Ink) is the 



eigenstate of the spectral number operator = a\dk- 

To illustrate this unraveling we consider two cases of a 
driven TLA coupled to the bath. In the first the 'bath' 
is a single mode harmonic oscillator and in the second 
the bath consists of three modes as described in the in- 
troduction. 



A. A single-mode bath 

The first example consists of a bath with only one mode 
{k — 1). Thus it only has one bath hidden variable with 
value n at time t. Thus n(t) denotes the trajectory of 



this value through time. By Eq. (2.10) the interaction 
Hamiltonian in the interaction frame with the free dy- 
namics removed is 



Vintit)=ih[g*e 



-i(wi-wo)t;it 



a' a\ 



(4.2) 



Choosing the bath frequency, lui = luq and including the 
driving of the TLA the evolution for the guiding state is 
determined by 



* t 

— (Tj; + g aa^ 



\^it)). (4.3) 



To calculate the stochastic evolution of the property 
state, \^n{t){t)) = \n{t))\tpn{t){'t)) , we use Bell's solution 
to find Jnm{t) [see Eq. (3.13)]. Doing this we get 

Jnra{t) = 2lu\{{-ipn{t)\{n\[ig* CFo) - zgo-t a] |m) | (t )) } . 

(4.4) 

where |?/;„(t)) = (?i|5'(t)) is the unnormalised conditioned 
state in configuration space. This simplifies to 



Jn,n{t) - 2Im[i5 (7/.„ (t ) I CT Vrn (0 > V?l + 1 (5 

Tl— 1,771 ]. (4.5) 

Here we see that for a given m only J^^i „i(t) arc non- 
zero. Thus the only transitions allowed arc from rn —t 
TO± 1. The transition matrices [see Eqs. (3.8) - (3.11)] 
are as follows. 

For J„i+i,,„(i) > (or J,„-i,m(t) < 0), 



T, 



7n+ 1,771 



J, 



7n+ 1,771 



(i)/Pr(z„,i) 



= 2Im[zg*(V'™+i(t)|<T|V™(t)) 

/(m + l)Pr(m+ l,t) 
Pr(m,i) 

2m-l,m(i) = 0, 

and for Jm+i,m(i) < (or J„i-i,m(i) > 0), 
0, 



^771 + 1, Tn(^) 
^rn — l,m(^) 



(4.6) 
(4.7) 

(4.8) 



~ <An— l,7n (^) /Pl"(-2^m 5 ^) 

= -2Im ig{il;„i-i{i)\cr^\i^yn{t)) 



I ?7iPr(m — l,t) 



(4.9) 



Using a random number generator we can simulate a 
typical trajectory for n{t) based on the above transition 
rates. For example given that n{t) is m at time t, depend- 
ing on which one of Jm±i.m(i) is positive [Jm+i,m{t) or 
Jm-i,m(t)] determines whether the value of the bath hid- 
den variable is going to jump up or down in the interval 
dt. Let us assume Jm+i,Tn{t) is positive then the proba- 
bility of an upward jump in the interval dt is Tm+i,m{t)dt, 
the simulation works by choosing a random variable be- 
tween [0, 1] and if Tm+i,m{i)dt is greater than this an 
upward jump occurs. This procedure is then repeated 
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until the desired simulation time is reached. Note in 
Eqs. (4.6) - (4.9) how the conditioned state in config- 
uration space, |'0„(t)), guides n{t). However, unlike dif- 
fusive non-Markovian we can see here, explicitly, that it 
is impossible to define this transition in terms of only one 
conditioned state, thus no analytical evolution equation 
for \ipn{t)(t)) can be derived. 

Choosing Q = 5g a numerical simulation for a typical 
trajectory is shown in Fig. 1 in Bloch representation 



x{t) = (&x)t, 

z{t) = {a:,)t, 



(4.10) 
(4.11) 
(4.12) 



where {A)t = (V'n(t)(i)|^|'0n(t)(^))- Also shown in this 
figure is the trajectory for the value of the hidden value, 
n{t). Here it is observed that the conditioned state 
evolves smoothly until there is a jump in n(t). Most 
of jumps are up in photon number (all except one which 
occurs at large t). This is because with classical driving 
we are effectively adding energy to the system and for 
it to dissipate this energy photons are emitted from the 
atom into the bath. This is further verified by the fact 
that for most jumps there is a corresponding lowering 
of atomic energy [that is, a lowering of z{t)]. Note for 
Markovian dynamics an upward jump in photon number 
always puts the system in the ground state (z = —1). 
However in this figure this is not observed and in fact 
(for large time in the figure) there are upward jumps 
in atomic energy when the value of n{t) increases. This 
clearly shows that non-Markovian dynamics is a lot more 
complicated than Markovian. Note the above strange be- 
haviour (not Markovian-like) becomes more evident at 
later times, when there are many photons in the bath. 
This is expected since in the Markovian case there is 
never more then one photon per bath mode on average. 

To show that by performing an ensemble average over 
these trajectories do in fact give the reduced state, the 
difference between the ensemble average of 1000 trajec- 
tories and the actual reduced state is shown in Fig. 2. 
Here it is observed that to within statistical error these 
methods agree. 



B. The three mode bath 
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FIG. 1: Tliis figure depicts a typical trajectory in Bloch rep- 
resentation for the spectral mode unraveling when the system 
is a driven TLA and the bath is a single mode harmonic os- 
cillator. In this figure all calculations were done using the ini- 
tial system state \tp{0)} = \b) with system parameters g = 1, 
n = 5. Time is measured in units g^^ . 
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time, t [g"''] 

FIG. 2: This figure depicts the diff'erence between the ensem- 
ble average of 1000 trajectories and the actual reduced state 
for the spectral mode unraveling when the system is a driven 
TLA and the bath is a single mode harmonic oscillator. All 
parameters are defined in Fig. 1. 



In the above example we observed that it is possible 
to derive numerically the equivalent to a non-Markovian 
SSE for the spectral mode unraveling for a simple system. 
Now we consider a slightly more complicated system: a 
TLA coupled linearly to a three mode optical bath, with 
a central frequency equal to the frequency of the atom, 
Wo, and the two side bands have frequencies uJozLQ. That 
is, the interaction Hamiltonian [Eq. (2.10)] becomes 



With this interaction Hamiltonian the guiding state evo- 
lution is determined by 



}l*(0>- 



(4.14) 



Vint(i) 



k=-l 



ikQ.t ~ 1 



(4.13) As in the last example, to calculate the stochas- 
tic evolution of the property state [|^'{„^.(t)}. (t)) = 



\{nk{t)})\ip{n^y(t))], we use Bell's solution to find 1 



J, 



n_ 1 ,m_-i ;no ,mo ;ni ,m 



[see Eq. (3. 13)]. Here we find that 
only six values of Jn_i,m_i;no,mo;ni,mi (0 arc non-zero. 

We expect that in the short time and large driving limit 
the conditioned state will jump between approximate ax 
eigenstates {\(Tx = 1) and Icr^: = —I))- This expectation 
can be illustrated analytically by considering the follow- 
ing argument. We move to a second interaction frame, 
this being the frame with both the free dynamics (stan- 
dard interaction frame) and i?int removed. That is, the 
interaction Hamiltonian is 



1 

Vint'(i) = X] {■9fc['^2; + 



k=~l 
-iOtl - 1 ^ikitt 



+a e lofee 



gkWx + o-Je ' 
}/2, (4.15) 



where 



at = |a, = l)(a, = -l|, 

K = \(^x = -i){<yx = i|, 



(4.16) 
(4.17) 



and the prime denotes rotation to a second interaction 
frame. 

We can now make the assumption that that for large 
driving all sinusoidally varying terms can be neglected 
(this is effectively a second RWA). Thus Eq. (4.15) be- 
comes 

Vint'if) = i7i|.gli(7^aLi -.g-iCTjTa-i +30(7^4 
-goaxao + gla^a\ -5i(T+ai|/2. 

(4.18) 

Moving back to the original interaction frame results in 
Vint(0 = i?i|gliO-+aLi£ 



g-ia^ a_ie 



iflt 



-gia^aie 



(4.19) 



That is, in the large driving limit Eq. (4.13) approxi- 
mates to the above Hamiltonian. Now since a~ and 
(7+ are lowering and raising operators for the ax basis 
we can conclude that when there is a transition from 
the \ax ~ 1) to \ax = —1) eigenstates there is a corre- 
sponding upward jump in the photon number for the 
higher-frequency mode or a lowering in the photon num- 
ber for the lower-frequency mode. The transition from 
the \ax = — 1) to \ax = 1) eigenstatc result in jumps of 
opposite nature. A jump in the photon number for the 
central-mode should not change the conditioned state. 

These predictions are verified in Figs. 3 (A) and (B) 
where we show simulations of the exact Bell dynamics 
with interaction Hamiltonian Eq. (4.13). Note the slight 
variation from the above prediction which we believe is 




(B) 
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FIG. 3: This figure depicts a typical trajectory in Bloch rep- 
resentation for the spectral-mode unraveling when the system 
is a driven TLA and the bath consists of 3 modes. In this fig- 
ure all calculations were done using the initial system state 
ji/'(0)) = (je) -I- j&))/\/2 with parameters g-i = go = 5i = 1, 
and — 20g. Time is measured in units g^^ . Part (A) is for 
< t < 10 and part (B) for W<t < 20. 



due to the fact that the driving is not infinite. The first 
figure shows a typical trajectory for the conditioned state 
for parameters g-i = go = gi = 1, and f2 = 20g for the 
first 10 units of time whereas the second figure shows the 
next 10 units of time. It is in this second period of time 
where downward jumps in the outer modes are observed. 

The prime motivation behind considering this 3-mode 
system is that in the strong driving limit the jump dy- 
namics should exhibit similar features to those observed 
in a spcctral-modc unraveling of a Markovian open quan- 
tum system. This is because in the large driving limit 
the fluorescence of a driving TLA takes the form of a 
MoUow spectrum [25, 33] [Fig. 4 shows this fluorescence 
spectrum (line graph) for Q equal to 20 times the spon- 
taneous emission rate, 7]. In the large driving limit the 
MoUow spectrum approaches becoming delta functions 
centered on the frequencies ujq — Q, loq, and cjQ + fl. Thus 
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CO- (Op [y] 

FIG. 4: This figure depicts the fluorescence spectrum for a 
Markovian driven TLA with i} equal to 20 times the sponta- 
neous emission rate, 7. 



one would expect that in this limit an unraveling in terms 
of spectral modes would have similar dynamics in both 
a non-Mar kovian and Markovian open quantum system. 
But what is the spectral-mode unraveling for a Marko- 
vian open quantum system. In Ref. [34] an unraveling 
of a driven TLA with a Markovian system-bath inter- 
action using filter cavities and photdetectors was intro- 
duced. We call this the Wiseman and Toombes filtered 
Markovian (WTFM) unraveling. (Note this unraveling is 
an example of an unraveling of a Markovian open quan- 
tum system which has non-Markovian conditioned state 
evolution for the atom) . In this unraveling the field emit- 
ted by the atom is detected using two cavities tuned to 
the sidebands of the Mollow triplet (these act as mode 
filters). Thus in the strong driving limit the WTFM un- 
raveling effectively corresponds to a spectral-three-modc 
unraveling. In Ref. [34] it was observed that in this 
limit the conditioned state for the TLA jumps between 
almost ax eigenstates. Note, the non-Markovian case dif- 
fers from the WTFM in that it is possible to get down- 
ward jumps in photon number. In the Markovian case 
this is not possible. However, as before, these anomalous 
events become more frequent as the modes become more 
populated. 



V. TEMPORAL-MODE UNRAVELING 

In the Markovian limit the simplest Markovian SSE is 
direct detection. This detection scheme involves mea- 
surements performed into temporal modes of the sys- 
tem (rather than the frequency modes) [13]. A temporal 
mode is best defined by considering the electromagnetic 
field operator E(r). The standard definition of this (ne- 
glecting polarization) is [36] 



E(r) 



k 



k 



(5.1) 



where Ufe(r) if the mode function for the k^^ frequency 
(spectral) mode. This can be rewritten as 



E(r) = -z^6tv;(r)+z^S,v,(r), 



(5.2) 



where v,- (r) labels a new type of mode with annihilation 
and creation operators br and b^. These new modes can 
be related to the frequency modes by 



aklr,k, 



(5.3) 



where j^.k are the elements of a unitary matrix 
(Et 7r,fc7r,A;' = 4,i;' )• Choosing jr,k to givc a discrete 
fourier transform, 

1 " 

br = dfe exp(— i27rTfc/K), (5.4) 

^k 

1 " 

Qk = hr cxp{i2iTTk/ k), (5.5) 

results in Vi- (r) having the functional form of a temporal 
mode. 

Since we have now defined an annihilation operator for 
the temporal mode it is possible to define the observable 
to which we wish to give an objective reality. This is the 
temporal-mode- number operator. It = b].br. Thus the 
preferred projective measure for this unraveling is 



{7r4 = {i,y,® 1{/„J)({/„J1}, 
where 1{/„J) = \In,) ® ... ® ]/«,) ® ... ® ]/«, 



(5.6) 



and \In^) 

is an cigcnstate of Ir- This set of operators {Ir} are the 
hidden variables for this unraveling. Each operator It 
takes one of its possible integer values {In^}- 

To provide a clearer picture of this unraveling we con- 
sider briefly the Markovian case. In the Markovian limit 
the number of modes becomes continuous (k —^ 00) and 
the system-bath coupling becomes flat {gk 1/7/2 ). 
As such we must define continuous (in frequency) anni- 
hilation operators a{uj) where lu ^ uj ^ ujo with uj being 
a continues variable labelling the frequency of the bath. 
This results in the temporal modes, denoted by fo(r), be- 
coming continuous in time and being related to a(iu) by 
the fourier transform 



b{r) = 
a{oj) = 



1 



/27r 
1 



a{uj) exp{—iQT)duj., (5.7) 
b{T) cxp{iujT)dT (5.8) 



Furthermore the interaction Hamiltonian, Eq. (2.10) 
(this is in the interaction frame with the free system and 
bath dynamics removed) becomes 



Vint{t)=in 



(a))]da). 
(5.9) 
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Assuming that ojo is large (this is vahd for optical situ- 
ations), with little error the lower limits of the integrals 
can be taken as oo [35]. This then allows us to rewrite 
Eq. (5.9) as 

Vint(i) = / [&6{t - t)b\T) - a^S{t - T)b{T)]dT. 



(A) 



(5.10) 

Here we see that at time t only one temporal mode is 
coupled to the system. This is precisely why the bath 
in a Markovian open quantum system only effects the 
system at one time (no memory effects). 

For a non-Markovian open quantum system the above 
delta correlations will not exists and as such the bath 
will have a memory. Writing Eq. (2.10) in terms of the 
discrete temporal modes [Eq. (5.4)] gives 

v,^t{t) = ihJ2Kit)»H ~ Cr{t)a^K], (5.11) 



where 



Crit) 



1 



9ke 



'iilk t-\-i2nTk/ K. 



(5.12) 



This clearly does not have a delta correlation between 
t and T as occurs in the Markovian case [Eq. (5.10)]. 
To show that this does occur in the Markovian limit we 
simply let gk become fiat (gk = g) and assume a constant 
spacing in frequency Qk = kfl. Doing this results in 

(k-1)/2 

c^(t) = -^(l + 2 V cos\k(nt-2TrT/K)]\. (5.13) 
t{ ^ 

This by definition is a Kroneckcr (5-function in the k 
oo limit. This in turn implies a Dirac (5-function in the 
continuous limit. 

To illustrate this unraveling we consider a driven TLA 
coupled to a three mode bath as before. For the three 
mode case Eq. (5.4) becomes 



1 ^ 

br = —!=■ dfc exp(— i27rTfc/3). 
V 3 



(5.14) 



Using these operators we can rewrite the interaction 
Hamiltonian, Eq. (5.11), as 

1 

Mnt(t) = ih Y,{'''r{t)abl~Cr{t)&^br]. (5.15) 



where 



Crit) 



9 

x/3 



|l + 2cos(m-27rT/3)|. (5.16) 



Here we have assumed g-i = go = gi = g as before. 
Although these functions are not ^-functions, c^(i) are 
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FIG. 5: This figure depicts a typical trajectory in Bloch rep- 
resentation for the temporal-mode unraveling when the sys- 
tem is a driven TLA and the bath consists of 3 modes. In 
this figure all calculations were done using the initial system 
state |'i/'(0)) = |&) with parameters g-i — go = gi = 1, and 
= 20g. Time is measured in units g~^. Part (A) is for 
<t <5 and part (B) for 5 < t < 10. 



peaked at times t ~ 2tt{t/3 + n)/Cl for n an integer, as 
shown in the finial plot of Fig. 5. 

Using the same techniques as in the last unravel- 
ing it is possible to determined a probability current 
Jni,mi-n2,m2;n3,m3it). Then we cau find a typical tra- 
jectory for the values of the bath hidden variables 
{In^{t)} and a numerical trajectory for the condition 
state \ip{i„^{t)}{t))- Figures 5 (A) and (B) illustrates 
a typical trajectory for a driven TLA with driving fre- 
quency il = 20g. Here we see that the dynamics are 
clearly unlike that observed for a Markovian open quan- 
tum system under a temporal-mode unraveling (direct 
detection). In the Markovian case an upward jump in 
the temporal mode corresponds to the conditioned sys- 
tem state jumping to the ground state {z = —1); see Ref. 
[13]. However for an upward jump into the r"^ temporal 
mode which occurs when c^(t) is a maximum this be- 
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haviour is observed, and the closer c^(t) is to maximum 
the closer the conditioned state jumps to z = —1. Thus 
we expect that when there are enough modes such that 
at every time t one of the c^{t)' s is a maximum (Marko- 
vian limit) an upward jump in photon number will re- 
sult in the system state jumping to the ground state. 
Note, as in all the other examples at later times when 
the modes become populated extra non-Markovian fea- 
tures are present. Firstly downward jumps in temporal 
mode number are present. These result in the condi- 
tioned state jumping towards the excited state (z = 1). 
But more interesting at approximately t = 4.9,9^^ (as 
well as at i = 9.5f;~^), there is a upward jump in the 
r = — 1 (r = 0) temporal mode which results in an in- 
creasing of atomic energy (the conditioned state jumping 
towards a higher atomic energy state) . 



VI. COMPARISON WITH JACK, COLLET, 
AND WALLS 

Recently Jack, Collet and Walls (JCW) have published 
a series of papers, Refs. [26, 27, 28], where they devel- 
oped a quantum trajectory theory (continuous in time 
measurement theory) for non-Markovian systems which 
incorporates jump-like unravelings. Thus the aim of this 
section is to review their theory and discuss how it differs 
from ours. 

In the JCW theory, because it concerns measurements, 
they have to add the requirement, "a measurement at a 
certain time does not disturb any future measurements" 
[28] . This basically means that the average state of the 
system (average over all possible measurements records) 
must be equal to the reduced state. Thus there is no 
back action on the system from the measurements. How 
is this possible? For argument's sake, consider two con- 
secutive quantum non-demolition measurements of the 
bath observable 

Z = {(^Zn, Isys ® \Zn) {Zji I)}. (6.1) 

Then we can write the average state of the system as 

Pit) = 5](4^)|J7(t2,ti)|zW)(zW|?7(ti,0)|*(0)) 

x(vI/(0)|J7t(ti,0)|zW)(4i)|,7t(i^^i^)|,(2)) 

(6.2) 

where the superscripts (1) and (2) refer to the first mea- 
surement made at time ii and the second measurement 
made at time ^2- This is only equal to the reduced state, 
Eq. (1.2), if [U{t2,ti), |z^^)(z^'|] = 0. This implies that 
either U{t2,ti) is diagonal in the bath basis or 

U{t2,ti) and {|zTO')(zm^|} operate on separate Hilbert 
spaces. This is precisely what happens in the Markovian 
limit (because of the delta correlations between the bath 
noise operators [35]). JCW theory uses the later case also 



and proposes that we can define a field quantity (j){t) to 
be measured, such that 

m,$\t')]=S{t-t'). (6.3) 

However, apart from Markovian open quantum systems 
it is not clear what measurement schemes this is appli- 
cable too. JCW consider in Ref. [26] a non-Markovian 
unraveling of a Markovian open quantum system, this 
being their spectral detection unraveling. This is similar 
to the WTFM unraveling, except they do not extend the 
system to include cavities. 

In the JCW theory to account for the the memory 
effects of non-Markovian unravelings it is imposable to 
assign a pure state for the system at time t conditioned 
on a past measurement record, l[o,t] (a string of results 
n, ...,rfc where t = kSt). Instead the state of the system 
conditioned on the measurement record l[o.t] is given by 

where 

Pr(I[o,t])= |d{^^>^}(^i,,,,,(<)l^i[o,.iW) (6.5) 

and 

Pi[o,,W= / d{r^>^}|^i[„„,(i)>(^i,o,.,*)l- (6-6) 

Here |V'i[o ,] (i)) is the the state of the system at time t 
conditioned on the complete record for all time (measure- 
ment which have not yet come to be). In a sense this is 
like a retrodictivc state and may have some interpreta- 
tion under retrodictivc quantum mechanics [37] . Thus we 
see that the best we can do is to assign a mixed state to 
the system, this being a mixture over all possible future 
records. 

However, as JCW point out, if we can assume a finite- 
time memory function (a memory function which is zero 
for time less then t — tm in the past), then we can as- 
sign a pure state to the system at time t — tm condi- 
tioned on measurements up and until time t. That is, 
the state of the system conditioned on the record I[o,f] is 

\^i[o t] (t ^ tm)) (when appropriately normalized). 

To summarize, while their theory seems to be correct 
we are not sure of the applicability of it. In general non- 
Markovian systems will have back action effects if a mea- 
surement is performed, so the average state for consec- 
utive measurements will not be the reduced state. Also 
the pure states they associate for the system arc not the 
same as our conditioned states. In fact in Ref. [8] we 
came to the conclusion that the only possible interpre- 
tation of non-Markovian SSEs (our conditioned states) 
under the orthodox theory is they are numerically tools 
used to generate the correct state of the system at time 
t given that a measurement has been performed on the 
bath at this time yielding the appropriate result. 



11 



VII. DISCUSSION AND CONCLUSION 

In this paper we have investigated non-Markovian 
unravehngs that exhibit jump-Hke behaviour. We ob- 
served that it is impossible to define an analytical non- 
Markovian SSE for these unravelings, but by using the 
modal interpretation of quantum mechanics, namely 
Bell's beable dynamics [22], it is possible to calculate the 
numerical equivalent to the solution of a non-Markovian 
SSE. To illustrate this we considered an open quantum 
system consisting of a TLA coupled to a bath of harmonic 
oscillators. 

The first example we considered was a driven two level 
atom (TLA), driven with Rabi frequency fl, and coupled 
linearly to a single mode bath. Here we observed that the 
non-Markovian dynamics (conditioned states evolution) 
are quite different in nature to Markovian dynamics. For 
example we observe that it possible to get upward jumps 
in photon number which result in an increasing of the 
atomic energy of the atom. 

The second example was again a driven TLA, but we 
considered the bath to contain three modes. The cen- 
tral mode of the bath had a frequency equal to the fre- 
quency of the atom, luq, while the two outer modes had 
frequencies loq ± O. This was chosen as the spectrum of 
this system (three modes) is similar to the fluorescence 
spectrum of a driven TLA (MoUow spectrum [25] ) in the 
strong driving and Markovian limit (see Fig. 4). Thus 
one expects that the conditioned states evolution should 
contain some features which are consistent with the dy- 
namics of a Markovian open quantum system. To illus- 
trate this we considered the two jump-like unravelings: 
the spectral mode and temporal-mode. It was observed 
that in both these unravelings the jump-like dynamics 
did contain some features of the equivalent Markovian 



open quantum system. 

In conclusion, by accepting the modal interpretation 
of non-Markovian SSEs it is possible to generalize non- 
Markovian unravelings to include jump-like unravelings, 
but it is only possible to numerically determine the so- 
lution. Note this generalization also applies to diffusive 
non-Markovian SSEs. That is, we can extend the co- 
herent [1, 2, 3, 8, 12], position [11, 12] and quadrature 
[7, 8, 11] unraveling to include all possible choices of a 
continuous preferred projective (and positive operator) 
measures. However, in the diffusive case it maybe pos- 
sible to derive a generalized analytical non-Markovian 
SSE. Bassi [10] has already proceeded down this path by 
calculating a generalized linear non-Markovian SSEs, for 
diffusive unravelings, but further question still remain. 
For example, what subset of reduced states do these lin- 
ear non-Markovian SSEs belong to and how can an ex- 
tension to the normalized case be made? Furthermore 
although it maybe possible to write an analytical expres- 
sion for a diffusive non-Markovian SSE we believe that, 
in general, due to the functional derivative the solution 
of this equation will only be determined by numerical 
perturbative techniques (see [4] or [9] for perturbative 
approximations) which in some circumstances effectively 
amounts to solving the Schrodinger equation for the total 
state l^'(i)). 
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